CATALYTIC BRANCHING PROCESSES VIA SPINE TECHNIQUES AND 

RENEWAL THEORY 



LEIF DORING AND MATTHEW ROBERTS 



Abstract. In this article we contribute to the moment analysis of branching processes in cat- 
alytic media. The many-to-few lemma based on the spine technique is used to derive a system 
of (discrete space) partial differential equations for the number of particles in a variation of 
constants formulation. The long-time behavior is then deduced from renewal theorems and 
induction. 



1. Introduction and Results 

A classical subject of probability theory is the analysis of branching processes in discrete or 
continuous time, going back to the study of extinction of family names by Francis Galton. There 
have been many contributions to the area since, and we present here an application of a recent 
development in the probabilistic theory. We identify qualitatively different regimes for the longtime 
behaviour for moments of sizes of populations in a simple model of a branching Markov process 
in a catalytic environment. 

To give some background for the branching mechanism, we recall the discrete-time Galton- Watson 
process. Given a random variable X with law /i taking values in N, the branching mechanism is 
modelled as follows: for a deterministic or random initial number Z € N of particles, one defines 
for n = 1,2, ... 

z n 

Z n +i = y^X r (n), 

where all X r (n) are independent and distributed according to fi. Each particle in generation n 
is thought of as giving birth to a random number of particles according to /z, and these parti- 
cles together form generation n + 1. For the continuous-time analogue each particle carries an 
independent exponential clock of rate 1 and performs its breeding event once its clock rings. 

It is well-known that a crucial quantity appearing in the analysis is to = E[X], the expected 
number of offspring particles. The process has positive chance of long-term survival if and only if 
to > 1. This is known as the supercritical case. The cases to = 1 (critical) and to < 1 (subcritical) 
also show qualitatively different behaviour in the rate at which the probability of survival decays. 
As this paper deals with the moment analysis of a spatial relative to this system, we mention the 
classical trichotomy for the moment asymptotics of Galton- Watson processes: 

(1.1) lim e-^-^ElZ?] e (0,oo) Vfc e N if to < 1, 

t— >oo 

(1.2) lim t^ElZ*] e (0,oo) Vfc e N if to = 1 

t— >oo 

(1.3) lim e-^-^ElZf] e (0,oo) Vfc e N if to > 1 

so that all moments increase exponentially to infinity if to > 1, increase polynomially if to = 1, 
and decay exponentially fast to zero if to < 1. 
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In the present article we are interested in a simple spatial version of the Galton- Watson process 
for which a system of branching particles moves in space and particles branch only in the presence 
of a catalyst. More precisely, we start a particle £ which moves on some countable set S according 
to a continuous-time Markov process with Q-matrix A. This particle carries an exponential clock 
of rate 1 that only ticks if £ is at the same site as the catalyst, which we assume sits at some fixed 
site € S. If and when the clock rings, then the particle dies and is replaced in its position by a 
random number of offspring. This number is distributed according to some offspring distribution 
and all newly born particles behave as independent copies of their parent: they move on S 
according to A and branch after an exponential rate 1 amount of time spent at 0. 

In recent years several authors have studied such branching systems. Often the first quantities 
that are analyzed are moments of the form 

M k (t,x,y)=E[N t (y) k \t; Q = x] and M k (t, x) = E[N k | £o = x] , 

where N t (y) is the number of particles alive at site y at time t, and N t — J2 y es Nt(y) is the total 
number of particles alive at time t. Under the additional assumption that A — A is the discrete 
Laplacian on Z d , the moment analysis was first carried out in [ABY98J, ABY98b , AB00 via 
partial differential equations and Tauberian theorems. More recently, the moment analysis, and 
moreover the study of conditional limit theorems, was pushed forward to more general spatial 
movement A assuming 

(Al) deducibility, 

(A2) spatial homogeneity, 

(A3) symmetry, 

(A4) finite variance of jump sizes. 

Techniques such as Bellman-Harris branching processes (see |TV03j . [VT05] . |B10j . |B11| ). operator 
theory (see |Y10j ) and renewal theory (see |HVT10p have been applied successfully. Some of these 
tools also apply in a non-symmetric framework. We present a purely stochastic approach avoiding 
the assumptions (A1)-(A4). In order to avoid many pathological special cases we only assume 

(A) the motion governed by A is irreducible. 

This assumption is not necessary, and the interested reader may easily reconstruct the additional 
cases from our proofs. 

In order to analyze the moments M k one can proceed in two steps. First, a set of partial differential 
equations for M k is derived. This can be done for instance as in [AB00| via analytic arguments from 
partial differential equations for the generating functions E x [e~ zNt ^] and E^le -2 ^*] combined 
with Faa di Bruno's formula of differentiation. The asymptotic properties of solutions to those 
differential equations are then analyzed in a second step where more information on the transition 
probabilities corresponding to A implies more precise results on the asymptotics for M k . This is 
where the finite variance assumption is used via the local central limit theorem. 

The approach presented in this article is based on the combinatorial spine representation of [ HRllj 
to derive sets of partial differential equations, in variation of constants form, for the fcth moments 
of N t (y) and Nt. A set of combinatorial factors can be given a direct probabilistic explanation, 
whereas the same factors appear otherwise from Faa di Bruno's formula. Those equations are 
then analyzed via renewal theorems. We have to emphasize that under the assumption (A) only, 
general precise asymptotic results are of course not possible so that we aim at giving a qualitative 
description. Compared to the fine results in the presence of a local central limit theorem (such as 
Lemma 3.1 of [HVTlOj for finite variance transitions on Z 4 ) our qualitative description is rather 
poor. On the other hand, the generality of our results allows for some interesting applications. 
For example, one can easily deduce asymptotics for moments of the number of particles when 
the catalyst is not fixed at zero, but rather follows some Markov process of its own, simply by 
considering the difference walk. 
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To state our main result, we denote the transition probabilities of A by Pt(x, y) = P x (£t — y) and 
the Green function by 



/>oo 

Goo{x,y) = / p t (x,y)dt. 
Jo 



Recall that, by irreducibility, the Green function is finite for all x, y £ S if and only if A is 
transient. For the statement of the result let us further denote by 



Lt{y) = [ l {6e = y} ds 
Jo 

the time of £ spent at site y up to time t. 

Theorem 1. Suppose that fj, has finite moments of all orders; then the following regimes occur 
for all integers k > 1: 

i) If the branching mechanism is subcritical, then 

lim M (t,x) G (0,oo) if A is transient, 

t— >oo 



t— »oo 

and 



lim M k (t,x) =0 if A is recurrent, 



lim M k (t,x,y) =0 in all cases. 



t— >oo 



ii) If the branching mechanism is critical, then 

M k (t x) 

lim / n \ fc _n € (0,oo) and M (t, x, y) = p t {x, y). 

iii) If the branching mechanism is supercritical, then there is a critical constant 

1 1>1 



" 600(0,0) 

such that 
a) for m < f3 

lim M 1 ^, x) € (0,oo) and lim M 1 ^, x, y) = 0; 

further, there exist constants c and C such that 

cE x [L t (0) k - 1 } < M k (t,x) < Ct k -\ 



b) for m = /3 



lim M k (t,x) = 00, 

t— >oo 



and 



/>oo 

lim M fe (i, a;, j/) = 00 if / rp r (0, 0) dr = 00 

M k (t x v) r°° 

lim e (0, 00) if / rp r (0, 0) dr < 00. 



(In both cases the growth is subexonential.) 
c) for m > /3 

lim e- kr< - m ^M k {t, x, y) e (0, 00) and lim e - kr ^ t M k (t, x) e (0, 00) 



where r(m) equals the unique solution A £0 J °° e A£ pt(0, 0) 



rn—l ' 
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We did not state all the asymptotics in cases ii) and iii)a). Our methods, see Lemma |3j do allow 
for investigation of these cases too; in particular they show how M k (t, x, y) can be expressed recur- 
sively by M l (t, x, y) for i < k. However, without further knowledge of the underlying motion, it is 
not possible to give any useful and general information. If more information on the tail of pt(x, y) 
is available then the recursive equations can indeed be analyzed: for instance for kernels on Z d 
with second moments the local central limit theorem can be applied leading to p t (x,y) ~ Ct~ d / 2 , 
and such cases have already been addressed by other authors. 

The formulation of the theorem docs not include the limiting constants. Indeed, the proofs give 
some of those (in an explicit form involving the transition probabilities pt) in the supercritical 
regime but they seem to be of little use. The use of spectral theory for symmetric Q-matrices A 
allows one to derive the exponential growth rate r(m) as the maximal eigenvalue of a Schrodinger 
operator with one-point potential and the appearing constants via the eigenfunctions. Our renewal 
theorem based proof gives the representation of r(m) as the inverse of the Laplace transform of 
Pt(0, 0) at l/(m — 1) and the eigenfunction expressed via integrals of pt(0, 0). As p*(0, 0) is rarely 
known explicitly, the integral form of the constants is not very useful (apart from the trivial case 
of Example 1 below). Only in case hi) b) for J Q rp r (0,0) dr = oo are the proofs unable to give 
strong asymptotics. This is caused by the use of an infinite-mean renewal theorem which only 
gives asymptotic bounds up to an unknown factor between 1 and 2. There is basically one example 
in which pt(0, 0) is trivially known: 

Example 1: For the trivial motion .4 = 0, i.e. branching particles are fixed at the same site as 
the catalyst, the supercritical cases iii) a) and b) do not occur as A is trivially recurrent so that 
13 = 1. Furthermore, in this example pt(0, 0) = 1 for all f > so that r(m) = m — 1. In fact by 
examining the proof of Theorem [l] one recovers ( 1.1|1.2|1~3 1 with all constants. 



The explicit representation for the exponential growth rate allows for a more careful comparison 
with the non-spatial case. 

Corollary 1. Let r(m) be the exponential growth rate obtained in the supercritical case of Theo- 
rem 1. Then 

m r(m) 

is convex, with 

rim) 

r(l) = 0, r(m) < m - 1, lim - v '— = 1. 

m->oo (m — 1) 

Proof. This follows from elementary manipulations of the defining equation for r(m). □ 

Remark 1. We reiterate here that our results can be generalized when the fixed branching source 
is replaced by a random branching source moving according to a random walk independent of the 
branching particles. For the proofs the branching particles only have to be replaced by branching 
particles relative to the branching source. 

2. Proofs 

The key tool in our proofs will be the many-to-few lemma proved in jHRllj which relies on modern 
spine techniques. These emerged from work of Kurtz, Lyons, Pemantle and Peres in the mid-1990s 
[KLPP971 IL971 ILPP95| . The idea is that to understand certain functionals of branching processes, 
it is enough to carefully study the behaviour of one special particle, the spine. In particular very 
general many-to-one lemmas emerged, allowing one to easily calculate expectations of sums over 
particles like 



E 



E /(«) 



IveNt 

where f(v) is some well-behaved functional of the behaviour of the particle v up to time t, and 
N t here is viewed as the set of particles alive at time t, rather than the number. It will always be 
clear from the context which meaning for N t is intended. 
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It is natural to ask whether similar results exist for higher moments of sums over N t . This is the 
idea behind ill HI lj. wherein it turns out that to understand the fcth moment one must consider 
a system of fc particles. The fc particles introduce complications compared to the single particle 
required for first moments, but this is still significantly simpler than controlling the behaviour of 
the potentially huge random number of particles in N t . 

While we do not need to understand the full spine setup here, we shall require some explanation. 
For each k > let — P(X — fc) and = E[A fc ], the fcth moment of the offspring distribution 
(in particular mi = m). We define a new measure Q = Qji, under which there are fc distinguished 
lines of descent known as spines. The construction of Q relies on a carefully chosen change of 
measure, but we do not need to understand the full construction and instead refer to |HR11| . In 
order to use the technique, we simply have to understand the dynamics of the system under Q. 
Under Q^i particles behave as follows: 

• We begin with one particle at position x which (as well as its position) has a mark fc. We 
think of a particle with mark j as carrying j spines. 

• Whenever a particle with mark j, j > 1, spends an (independent) exponential time with 
parameter rrij in the same position as the catalyst, it dies and is replaced by a random 
number of new particles with law Aj . 

• The probability of the event {Aj — a} is a?p a mj . (This is the jth size-biased distribution 
relative to /i.) 

• Given that a particles v%, . . . ,v a are born, the j spines each choose a particle to follow 
independently and uniformly at random. Thus particle Vi has mark I with probability 
(1 — a^ 1 )^ 1 , I = 0, . . . , j, i — 1, . . . , a. We also note that this means that there are 
always fc spines amongst the particles alive; equivalently the sum of the marks over all 
particles alive always equals fc. 

• Particles with mark are no longer of interest (in fact they behave just as under P, 
branching at rate 1 when in the same position as the catalyst and giving birth to numbers 
of particles with law [i, but we will not need to use this). 

For a particle v, we let X v (t) be its position at time t and B v be its mark (the number of spines it 
is carrying). Let a v be the time of its birth and t v the time of its death, and define a v (t) — a v At 
and T v (t) — t v A t. Let \t be the current position of the zth spine. We call the collection of 
particles that have carried at least one spine up to time t the skeleton at time t, and write skel(f). 
Figure [T] gives an impression of the skeleton at the start of the process. 



FIGURE 1. An impression of the start of the process: each particle in the skeleton 
is a different colour, and particles not in the skeleton are drawn in pale grey. The 
circles show the number of spines being carried by each particle in the skeleton. 

A much more general form of the following lemma was proved in [HRllj . 
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Lemma 1 (Many-to-few). Suppose that f 



is measurable. Then, for any k > 1, 



E 



E /(*«x(*)) 

...,v k eN t 



f(X Vh (t)) 



f(xl)---f(x k t ) II exp (m B „-l) 

i)Gskel(t) V 



«(*) 



,(t) 



l {X v (s))ds 



Clearly if we take / = 1, then the left hand side is simply the fcth moment of the number of 
particles alive at time t. The lemma is useful since the right-hand side depends on at most k 
particles at a time, rather than the arbitrarily large random number of particles on the left-hand 
side. 

Having introduced the spine technique, we can now proceed with the proof of Theorem 1. We first 
use Lemma[T]for the case k = 1, which is simply the many-to-one lemma, to deduce two convenient 
representations for the first moments: a Feynman-Kac expression and a variation of constants 
formula. Indeed, the exponential expression equally works for other random potentials and, hence, 
is well known for instance in the parabolic Anderson model literature. More interestingly, the 
variation of constants representation is most useful in the case of a one-point potential: it simplifies 
to a renewal type equation. Understanding when those are proper renewal equations replaces the 
spectral theoretic arguments of |ABY98j and explains the different cases appearing in Theorem 1. 

Lemma 2. The first moments can be expressed as 

(2.1) M 1 {t,x)=E x [e( m -Vti lo ^ dr ], 

(2.2) M l {t, x, y) = E, [e^ £ , 

where £ t is a single particle moving with Q-matrix A. Furthermore, these quantities fulfill 

(2.3) M 1 (i, x) = 1 + (m-l)p t {x,0) * M x (i, 0), 

(2.4) M\t, x, y) = p t (x, y) + (m- l)p t (x, 0) * M\t, 0, y), 
where * denotes ordinary convolution in t. 

For completeness we include a proof of these well-known relations. First let us briefly mention 
why the renewal type equations occur naturally. The Feynman-Kac representation can be proved 
in various ways; we derive it simply from the many-to-few lemma. The Feynman-Kac formula 
then leads naturally to solutions of discrete-space heat equations with one-point potential: 

J §- t u(t, x) = Au(t, x) + (m- l)l (x)u(t, x) 
\u(Q,x) = l y {x) 

Applying the variation of constants formula for solutions gives 



u(t,x) = Ptu(0,x) + / Pt- S (m — l)ta(x)u(s, x) ds 



Pt{x,v) + (m- 1) / p t - s (x,0)u(s,x)ds, 
Jo 



where P t is the semigroup corresponding to A, i.e. P t f(x) = E x [/(£ t )] 



Proof of Lemma^ To prove ( |2.1[ ) and (2.2) we apply the easiest case of Lemma [TJ we choose 
k = 1 and / = 1 (resp. f(z) = l y (z) for (2.2)). Since there is exactly one spine at all times, 
the skeleton reduces to a single line of descent. Hence m^,, — 1 = m — 1 and the integrals in the 
product combine to become a single integral along the path of the spine up to time t. Thus 

M\t,x) =Q x [e( m - 1 )/oio(«^ and M ^t,x,y)=Q x [e^foM^dr lyi ^ 
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which is what we claimed but with expectations taken under Q rather than the original measure 
P. However we note that the motion of the single spine is the same (it has Q-matrix A) under 
both P and Q, so we may simply replace Q with P, giving (2.1) and (2.2). 

The variation of constants formulas can now be derived from the Feynman-Kac formulas. We only 
prove the second identity as the first can be proved similarly. We use the exponential series to get 



= E 



L(m-l)/o lo 



E 



(m - 1)" 



E 



(m - l) n rt 



lo(fn)"-lo(fr n )<frn...dril v (&) 



Pt{x,y) +E 2 



OO r-t ft r-t 

J2(m-ir ■■■ l ^ ri )---l (^Jdr n ...dr 2 dr 1 l y ^ t ) 

n=l ''1 •'Tn—l 



The last step is justified by the fact that the function that is integrated is symmetric in all 
arguments and, thus, it suffices to integrate over a simplex. We can exchange sum and expectation 
and obtain that the last expression equals 

rt oo 

Pt (x,y) + (m-l) ^(m-l)"- 1 / .../ P K [£ ri = 0, . . . , & n = 0] dr n . . . dr 2 dn. 

Jo n=1 Jri Jr„_i 

Due to the Markov property, the last expression equals 



Pi 



(x,y) + (m-l) p ri (a;,0)^(m-l) n - 1 / .../ PoKm-n = 0, . . . , Zr n - ri = 0] dr n . . . dr 2 dn 

Jo n = i Jri Jr n -i 



and can be rewritten as 

r-l 



p t (x,y) + (m-l) / Pri (x,0) V(m-l)"- 1 
Jo V«_i 



t-n 



t-n 



, £r„ = 0] dr„ ...dr 2 dr\ . 



Using the same line of arguments backwards for the term in parentheses, the assertion follows. □ 

Having derived variation of constants formulas, there are different ways to analyze the asymptotics 
of the first moments. Assuming more regularity for the transition probablities, this can be done 
as sketched in the next remark. 



Remark 2. Taking Laplace transforms C in t, one can transform (2.3), and similarly (2.4-), into 
the algebraic equation 

£M\\,x) = - + (to — 1)£M 1 (A, 0)Cpx(x, 0) ,A > 0, 
A 

which can be solved explicitly to obtain 
(2.5) CM l (\x) 



, A > 0. 



A(1-(to-1)£ Pa (x,0)) 

Assuming the asymptotics ofp t (x,0) are known for t tending to infinity (and are sufficiently reg- 
ular), the asymptotics of Cp \{x, 0) for A tending to zero can be deduced from Tauberian theorems. 
Hence, from Equation (2.5) one can then deduce the asymptotics of CM 1 (X,x) as A tends to zero. 
This, using Tauberian theorems in the reverse direction, allows one to deduce the asymptotics of 
M 1 (t,a;) fort tending to infinity. 

Unfortunately, to make this approach work, ultimate monotonicity and asymptotics of the type 
P t(x,0) ~ Ct~ a are needed. This motivated the authors of [ABY98] to assume (A4) so that by 
the local central limit theorem 

d/2 



Pt(x,0) 



t 



-d/2 
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As we did not assume any regularity for p t , the aforementioned approach fails in general. We 
instead use an approach based on renewal theorems recently seen in |DS10| . 

Proof of Theorem^ j "or M 1 . Taking into account irreducibility and the Markov property of A, 
we see that the property "f^° lo(£ r ) dr — oo almost surely" does not depend on the starting 



value £o- To prove case i), we simply apply dominated convergence to (2.1) and (2.2). If A is 



transient, then l (£ r ) dr < oo almost surely and M 1 (t, x) converges to a constant. On the 
other hand if A is recurrent, then lo(£r) dr = oo almost surely and M x (t, x) — > 0. In both cases 
M 1 (i,x,y) — > 0, because if A is transient then l{j t=y } — > almost surely, and if A is recurrent 
then M 1 (t, x, y) < M x (t, x) 0. 

Regime ii) is trivial as here M 1 (t 1 x) ~ 1 and M 1 (t,x,y) = pt(x,y). Next, for regime hi) a) we 
exploit both the standard and the reverse Holder inequality for p > 1: 

(2.6) M^t, x, y) > E x [ e -(V(p-i))(Tr.-i) /<? M60 y )P, 

(2.7) M^t.z.y) <K x [e^ m - l ^oMir)dr^/p p ^ x ^( P -i)/ P _ 

In the recurrent case Go^O,!]) = 00 and thus f3 = 1, so this case has already been dealt with 
in regime ii). Hence we may assume that A is transient so tha t J n lo(£ r ) dr < 00 with positive 
probability. This shows that the expectation in the lower bound ( |2.6[ ) converges to a finite constant. 
By assumption m— 1 < /? so that there is p > 1 satisfying p(m — l) < p. With this choice of p, part 



3) of Theorem 1 of [DSlOj implies that also the expectation in the upper bound (2.7) converges 
to a finite constant. In total this shows that 

Cp t {x,yY < M\t,x,y) < C'pt^y)^ 1 ^ 

and the claim for M 1 (f, x, y) follows. For M 1 (i, x) we can directly refer to Theorem 1 of [DS10J. 

For regimes iii) b) and c) we give arguments based on renewal theorems. A closer look at the 



variation of constants formula (2.4) shows that only for x = 0, M 1 (t,a;,2/) occurs on both sides 
of the equation. Hence, we start with the case x = and afterwards deduce the asymptotics for 
x ^ 0. 

Let us begin with the simpler case iii) c). As mentioned above, in this case we may assume that A 
is transient so that J Q p r (0,0)dr < 00. Hence, dominated convergence ensures that the equation 
Jp 00 e~ xt pt(0, 0) dt = l/(m — 1) has a unique positive root A, which we call r(m). The definition 
of r(rn) shows that U(dt) := (m - l)e- r (" l )*p*(0, 0) dt is a probability measure on [0, 00) and 
furthermore e~ r ( mS>t pt(0,y) is directly Riemann integrable. Hence the classical renewal theorem 
(see page 349 of |F71j ) can be applied to the (complete) renewal equation 

f(t)=g(t) + f*U(t), 

with f(t) = e- r ( m ' t M 1 (t,0,?/) and g(t) = e~ r( - m ^p t (0,y). The renewal theorem implies that 

rf° q(s) ds 

(2-8) hm f(t) = J° y[ ' e (0, 00) 

*->oo J Q U((s, 00)) ds 

so that the claim for M 1 (t, 0,y) follows including the limiting constants. 

For iii) b), we need to be more careful as the criticality implies that (m — 1) p r (0, 0) dr = 1. 
Hence, the measure U as defined above is already a probability measure so that the variation of 
constants formula is indeed a proper renewal equation. The renewal measure U only has finite 
mean if additionally 



(2.9) / rp r (0,0)dr < 00 



In the case of finite mean the claim follows as above from (2.8| without the exponential correction 
(i.e. r(m) — 0). Note that pt(0,y) is directly Riemann integrable as the case j3 > implies that 
A is transient and Pt(Q,y) is decreasing. 
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If (2.9| fails, we need a renewal theorem for infinite mean variables. Iterating Equation (2.4 1 
reveals the representation 

(2.10) M\t, 0, y) = Pt (0, y)*J2( m ~ irMO, 0)*", 

n>0 

where *n denotes n-fold convolution in t and pt(0, y) *f>t(0, 0)*° = Pt(0, y). Note that convergence 
of the series is justified by 



(m-l)>(0,0)*" < l(m-l) J p r (0,0)dr 
and the assumption on m. Lemma 1 of E73J now implies that 



(2.11) 



]T(m-l)> t (0,0r 



n>0 



(m-l)Jl;J™Pr(0,Q)drds 



which tends to infinity as (m — 1) p r (0, 0) dr —> for s —> oo since we assumed that (m — 
l)p r (0,0) is a probability density in r. To derive from this observation the result for M 1 (t,Q,y) 1 
note that the simple bound pt(0,y) < 1 gives the upper bound 

(2.12! 

n>0 



M%0,y) < f £(m- 1)^(0, 0) w dr. 

JO „^ n 



For a lower bound, we use that due to irreducibility and continuity of p t (0,y) in t, there are 
< to < t\ and e > such that pt(0, y) > e for to < t < t\. This shows that 

rt—ti 

(2.13) 



M 1 (t,0,y)>e f ^(m-lTPrM^dr. 

■It-ta „>o 



Combined with (2.111 the lower and upper bounds directly prove the claim for M 1 (t, 0, y). 

It remains to deal with regime iii) b) and c) for x ^ 0. The results follow from the asymptotics 
of the convolutions as those do not vanish at infinity. But this can be deduced from simple upper 
and lower bounds similar to (2.121 and ( 2.13| ). 

The asymptotic results for the expected total number of particles M 1 (t,a;) follow from similar 
ideas: estimating as before 



M 1 (r,0)dr<M 1 (t,x)<l+ M 1 (r,0)dr 



t-t 



and applying case 2) of Theorem 1 of [DS10] to (2.1 1 with x — 0, the result follows. 



□ 



We now come to the crucial lemma of our paper. We use the the many-to-few lemma to reduce 
higher moments of N t and N t (y) to the first moment. More precisely, a system of equations is 
derived that can be solved inductively once the first moment is known. This particular useful 
form is caused by the one-point catalyst. A similar system can be derived in the same manner in 
the deterministic case if the one-point potential is replaced by a n-point potential. However the 
case of a random n-point potential is much more delicate as the sources are "attracted" to the 
particles, destroying any chance of a renewal theory approach. 

Lemma 3. For k > 2 the kth moments fulfill 

(2.14) M fc (t, x) = M\t, x) + M\t, x, 0) * g k ((M x (t, 0), • • • , M k ~\t, 0)) , 

(2.15) M fe (t, Xl y) = M\t 7 x, y) + M 1 (t, x, 0) * g k (M 1 (t, 0, y), ■ ■ ■ , M fe ~ 1 (t, 0, y)) , 
where 



g k (M\...,M k - 1 )=J2® 

J=2 



E 



i!,...,ij>0 
ii + ...+ij=k 



AP 1 ■ ■ ■ M lj . 
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Proof. We shall only prove equation (2.14); the proof of equation (2.15) is almost identical. We 



recall the spine setup and introduce some more notation. To begin with, all k spines are carried 
by the same particle £ which branches at rate m k = E[X k ] when at 0. Thus the k spines separate 
into two or more particles at rate m k — m when at (since it is possible that at a birth event all fc 
spines continue to follow the same particle, which happens at rate m). We consider what happens 
at this first "separation" time, and call it T. 



Let i\, . . . , ij > 0, i\ + . . . +ij = k, and define A k {j; i\, 



> 3 



to be the event that at a separation 



event, i\ spines follow one particle, ?2 follow another, . . . , and ij follow another. The first particle 
splits into a new particles with probability a^PaVa^ 1 (see the definition of Q fc ). Then given that 
the first particle splits into a new particles, the probability that i\ spines follow one particle, i^ 
follow another, . . . , and ij follow another is 

1 fa\ k\ 

n fe 



1 1 . 



.i 



(the first factor is the probability of each spine making a particular choice from the a available; 
the second is the number of ways of choosing the j particles to assign the spines to; and the third 
is the number of ways of rearranging the spines amongst those j particles). Thus the probability 
of the event A k (j; it, . . . , ij) under Q fc is 



1 



-E 



k\ 



.i ' 



(Note that, as expected, this means that the total rate at which a separation event occurs is 

k 



■ — Ve 



J=2 



E 



ii,...,ij>0 
ii + ...+ij=k 



i x \---ij\ 



since the double sum is just the expected number of ways of assigning k things to X boxes without 
assigning them all to the same box.) 

However, for j > 2, given that we have a separation event, Ak(j; ii, 



ij) occurs with probability 



1 



-E 



m k 



k\ 



m k 

i±l ■ ■ ■ ij\ \m k — m / 

Write Xt f° r the position of the particle carrying the k spines for t G [0, T), and define F t to be 
the filtration containing all information (including about the spines) up to time t. Recall that the 
skeleton skel(t) is the tree generated by particles containing at least one spine up to time t] let 
skel(s; t) similarly be the part of the skeleton falling between times s and t. Using the many-to-few 
lemma with / = 1, the fact that by definition before T all spines sit on the same particle and 
integrating out T, we obtain 



E[N t k ] 



TJ e (™s„-l) f£$ loCX„(s))<fa 
ueskel(t) 



e (m fc -l)/ T l (Xa)d.s 1{ 



7 ( , Yl e (m ^- 1)f ^ lo{XAs))ds 

i)£skel(T; t) 

+ Q k le^x- 1 ) lo(Xs)ds l {T>t} 

e (m k -l) /» lo(Xa)d S ( TOfe _ m )l ( Xu ) e -( m *- m ) fa lo(Xs)ds 



du 



t;Gskcl(ii; t) 

_j_ Q k e (TOfe-l)Jo £ lo(Xs)ds e -(m k -m) J* l (x 8 )ds 



CATALYTIC BRANCHING PROCESSES VIA SPINE TECHNIQUES AND RENEWAL THEORY 



11 



To prove equation (2.15), the same arguments are used with / = l y in place of / = 1. Now we 
split the sample space according to the distribution of the numbers of spines in the skeleton at 
time T. Since, given their positions and marks at time T, the particles in the skeleton behave 
independently, we may split the product up into j independent factors. Thus 



E[N*] = 



t k 

E E * 

j—2 ii,...,ij>0 
i 1 + ...+i j =k 

j 



i=i 



t k 

E E ^ 

j=2 i u ...,i s >0 
i 1 + ...+ij=k 



[(•): 


kl ink 




hi-- 


. ,M 




"[ n 

-u£skel(i — u 






+ 


Q k 






k 




E x [ 


[(■): 


hi-- 





,("»-!) / "lo(x.)d« 



MXu) 



(III 



,{m-l)f c ;i (Xs)d 3 



[jv„(o)] n ? " [^-j du + e - i N *\ - 



1=1 



where we have used the many-to-few lemma backwards with / = lo (first expectation) and / = 1 
(two last expectations) to obtain the last line. This is exactly the desired equation (2.14). For 
Equation (2.15) we again use / = t y in place of / = 1 and copy the same lines of arguments. □ 

Remark 3. The factors appearing in gk are derived combinatorially from splitting the spines. In 
Lemma 3. 1 of [ABOOj they appeared from Fad di Bruno 's differentiation formula. 

We need the following elementary lemma before we can complete our proof. 

Lemma 4. For any non-negative integer-valued random variable Y , and any integers a > b > I, 

E[y a ]E[Y"] > E[Y b }E[Y a - b+1 }. 

Proof. Assume without loss of generality that b > a/2. Note that for any two positive integers j 
and k, 

fk+jk a - j b k a - b+1 - j a - b+1 k b = jk(j - k) 2 (f- 3 + 2j a - 4 k + 3f- 5 k 2 + ... 

+ (a- b)j b - 2 k a - b - 1 + (a - b - l)j b - 3 k a - b + ... 
+ 2jk a - 4 + k a ^ 

> 0. 



Thus 



E[r a ]E[r] - E[r b ]E[r a - h+1 ] 

= j ap ( Y = j) E fcp ( y = fc ) ~ E ■ ? ' bp ( y = E ka ~ b+lp ( Y = k ) 

i>l k>l j>l 



k>l 



J2Y,U ak + k *i - / fca_b+1 - 3 a -" +1 k b W(Y = i)P(Y = k) 



j>l k>j 
> 

as required. 

We can now finish the proof of the main result. 



□ 



Proof of Theorem 1 for M k . Case i) follows just as for M 1 , applying dominated convergence to 
the Q fc -expectation in LemmaJT] Note that if T is the first split time of the k spines (as in Lemma|3]) 
then e( mh ~^ Jo 1 o(£s) ds [ s stochastically dominated by e ( m fc- 1 ) r where r is an exponential random 
variable of parameter — ra\\ this allows us to construct the required dominating random 
variable. 
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For case ii), using Lemmas [2] and [3] we find the lower bound 

(2.16) M k (t,x) > l + C [ p s (x,0)M k ~ 1 (t- s,0)ds >C I p s (x, 0)M k ~ 1 {t - s, 0)ds. 



An upper bound can be obtained by additionally using Lemma [4] (to reduce to the leading term 
M 1 M k - 1 ) to obtain 

(2.17) M k {t,x) < 1 + C [ Ps (x,0)M k - 1 (t- s,0)ds. 



Using inductively the lower bound (2.16) and furthermore the iteration 

v x (x Sl = o) / 1 ¥ Q (x S2 = o) . . . / 1 ^ 2 r (x Sk ^ = 0)ds fe _i . . . ds 2 d Sl 







(2.18) 



f f ... / P X (X S1 = 0, X S2 = 0, . . . , X Sk _, = 0)ds fc _i . . . ds 2 d Sl 

JO J Si J Sfc-2 

„ * , / V x (X sl =0,X 32 =0,...,X Sk _ 1 =0)ds k - 1 ...ds 2 dsi 
(k - 1)! Jo Jo Jo 



(k-iy. 
i 







t v * 

t{x 3 =o}ds 



we see that M k (t,x) goes to infinity if A is recurrent and to a constant if A is transient. This 
implies that the additional summand 1 in (2.171 can be omitted asymptotically in both cases. The 
claim follows. 

The lower bound of case iii)a) follows by the same argument as for case ii), and the upper bound 
is a straightforward induction using Lemmas [3] and |4j The cases ih)b) and c) also follow from 
Lemma [5] and induction based on the asymptotics for M 1 . □ 
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